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Abstract 

Dynamics and thermodynamics of the model with local anharmonism in the case of 
absence of the electron (Hubbard) correlation is investigated in the present work. The 
correlation functions, mean values of pseudospin and particle number as well as the ther- 
modynamical potential are calculated. The calculation is performed by diagrammatic 
method in the mean field approximation. Single-particle Green functions are taken in the 
Hubbard-I approximation. 

1 Introduction. 

The model considering the interaction of electrons with the local anharmonic mode of lattice 
vibrations is used in the last years in the theory of high-temperature superconducting crystals. 
Particularly, such property is characteristic for the vibrations of the so-called apex oxygen ions 
O/y along c-axis direction of the layered compounds of YBa2Cu307-type structure (see, [1-3]). 
An important role of the apex oxygen and its anharmonic vibrations in the phase transition into 
superconducting state has been already mentioned [4,5] and the possible connection between 
the superconductivity and lattice instability of ferroelectric type in high-T c superconducting 
compounds is under discussion [6,7]. In the case of local double-well potential, the vibrational 
degrees of freedom can be presented by pseudospin variables. The Hamiltonian of the derived 
in this way pseudospin-electron model has the following form [8] 

H = J2 H i + J2 tijKabj* , (1) 

i ijcr 

and includes besides the terms describing electron transfer (~ t^-), the electron correlation 
(£/-term), interaction with anharmonic mode (g-term), the energy of the tunnelling splitting 
(J7-term) and energy of the anharmonic potential asymmetry (/i-term) in the single-site part 

Hi = Un^nn + E (n^ + n tl ) + g(n^ + n^S? - SIS? - hS* . (2) 

Here, E Q gives the origin for energies of the electron states at lattice site (E Q = — //). 
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In this paper, which consists of two parts, our aim is to obtain the expressions for correlation 
functions which determine the dielectric susceptibility, mean values of pseudospin and particle 
number operators as well as the thermodynamical potential in the case Q = and absence of 
the Hubbard correlation U — 0. 

In the second part of the paper we shell perform numerical calculations for the analytical 
expressions obtained in the first part. We shall investigate values of pseudospin and particle 
number operators with the change of the asymmetry parameter h (T = const) or with the 
change of temperature T (h = const) for the cases of the fixed chemical potential value (regime 
\i = const) and constant mean value particle number. We shell present also some result for the 
case T = 0. On the basis of the obtained results the applicability of the Hubbard-I aproximate 
will be discussed. 



2 Hamiltonian and initial relations. 

We shall write the Hamiltonian of the model and the operators which correspond to physical 
quantities in the second quantized form using the operators of the electron creation (annihila- 
tion) at the site with the certain pseudospin orientation 



(3) 



6^(1/2 + , a+ =6+(l/2 + Sf) 
a oi = b oi (l/2-Sf) , a+ = 6+ (1/2 - Sf) 

Then we obtain the following expression for the initial Hamiltonian 

H = E{e("iT + n n) + e("tT + n n) - hS i}+ 

i 

+ EUj(a+a ja + at a a ja + ata ja + a+a ja ) = (4) 

ija 

— H Q + Hi n t, 

where 

e = E + g/2,e = E -g/2, (5) 

are energies of the single-site states; H Q is the single-site (diagonal) term, H int is the hopping 
terms. 

The introduced operators satisfy following commutation rules 

{a+, 5jv} = 8^8^(1/2 — Sf) , {a+,a jCT /} = 0, 
{a+ , a ja ,} = <^<W(l/2 + 5?) , {a+ , a ja ,} = . 

In order to calculate pseudospin mean values we shall use the standard representation of the 
statistical operator in form 

f f 

<j{(3) = T T exp < / H int (T)dr 



(8) 



Which gives the following expressions for (S, 



(s;) = -r^r^m)o = {s?mr ■ (9) 
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Here, the operators are given in the interaction representation 

A{t) = e TH °Ae- rH ° ; (10) 

the averaging (. . .) G is performed over the statistical distribution with Hamiltonian H , and 
symbol (. . .) c denotes the keeping of the connected diagrams. 



3 Perturbation theory for pseudospin mean values and 
diagram technique. 

The expansion of the exponent in (Q) in powers of H int (£D leads, after the substitution in 
equation (|^), to the expression that has the form of the sum of infinite series with the terms 
presented by the averages of the T-products of the electron creation (annihilation) operators 
at the site with the certain pseudospin orientation in the interaction representation. The 
evaluation of such averages can be performed using Wick's theorem. 

In our case this theorem has some differences from the standart formulation. Namelly, result 
of each pairing of operators (^) contains an operator factors, i.e. 



a 4 (r>+(r)= g(r> - r)5 io P^ 5 4 (r')5 + (r)= ~g(r> - r)5 io Pr, (11) 



at(T)a, t (T>)= -g(r' - r)5 io P t +, a+(r)5 i (r')= -~g{r' - r)5 io Pr. 

Finally, this gives the possibility to express result in term of the products of nonperturbated 
Green functions 

_ = (TMr)aiW)). = ^ ( ^ ; t > < 
{{aia+}) { tt^ ■ t >t, 

_ (T T a t (r)a+(rQ) _ i(T ,_ T) f :r>r', 

9io\T T ) - ,r~ ~ + - e °oi\ -l . / 

9io(r - t') = g(r - r')5 io , g io (r - t') = g(r - r')5 io , 
and averages of the certain number of the projection operators 

Pt = \ + S?, pr = l- S f. (13) 

Let us demonstrate this procedure for the case of evaluation of {Sf), for one of the terms which 
appear in the fourth order of the perturbation theory 

P P P P 

J dn J dr 2 J dr 3 J dn J2 ^AiiW^ais* ( 14 ) 
oooo ijhh -hhhh 

x (T T 5f af (ri)aj-(ri)aj {r 2 )a h (r 2 )a+ (r 3 )S i2 (r 3 )a+ (r 4 )a h (n)) . 
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The stepwise pairing of the certain operator with the other ones gives the possibility to 
reduce expression (0) to the sum of the averages of the smaller number of operators 

(r r 5fa+(r 1 )a i (r 1 )5+(r 2 )a il (r 2 )a+(r3)5 j2 (r 3 )a+(r4)a i 3(r4)) = 



(T T 3f at (n) a i (ri)a+ {T 2 )a h (r 2 )a+ (r 3 )a j2 (r 3 )a+ (r 4 )a i3 (r 4 ) ) + 



+ (T T Sf 0^^1)0^x1)5^^2)0^ (r 2 ) ai(r3)a i2 (r 3 )at(r4)a i3 (r4)> 



-<fo 3 ( r i - r 4) (T T 5f Pt aj -(ri)aJ (r 2 ) 0il (r 2 )a+ (r 3 )a i2 (r 3 )a+ (r 4 )) 



- ^)(r T 5fPta i (ri)at(r2)at(r 3 )a i2 (r3)at(r4)a i2 (r3)) . 



(15) 



The successive applications of the pairing procedure for ( |15"D leads, finally, to 

-9ih(n - r 2 )g ilj2 (r 2 - r 3 )ft 3j {r 4 - n) <7 i2jl (r 3 - r 2 ) ( T T S* P t P t P^r Pt ) - 

-^i3(n- r 4)^ 1 i 2 (r 2 -r 3 )^ 2 /r 3 -r 1 )^ 3il (r 4 -r 2 )(T r SfP+P+^ 
+^i 3 ( r i - TA)g ilj2 (T 2 - r 3 )g i2jl (r 3 - r 2 )g i:i] {r 4 - n) (T T 5f PtpJp-p+) . 

We introduce the diagrammatic notations 



(16) 



■| '\/\y\/\y\/^ 1 ' - £ ; 



-(gP^+gP;) o -s, z , 



and diagrams 




O 




o 




correspond to the expression ([H]). 
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The expansion of ([H]) in semi-invariants leads to the multiplication of diagrams (semi- 
invariants are represented by the ovals surrounding corresponding vertices with diagonal oper- 
ators and contain the 5-symbol on site indexes). For example, 




We shall neglect diagrams of types 2, 3, 5 i.e. the types including semi-invariants of the higher 
than first order in the loop (this means that chain fragments form the single-electron Green 
functions in the Hubbard-I approximation) and also the connection of two loops by more than 
one semi-invariant (this approximation means that selfconsistent field is taken into account in 
the zero approximation). 

Let us proceed to the momentum-frequency representation in the expressions for the Green 
functions determined on finite interval 0<t</3 when they can be expanded in Fourier series 
with discrete frequencies 

P n 

s(r) = ^£e^K); (17) 

P n 

1 . . 1 2n + 1 
; g(u n ) = -; u n = — - — vr . 

iuj n — e iuj n — e p 

The characteristic feature of the already presented diagrams and diagrams that correspond to 
other orders of the perturbation theory is the presence of chain fragments. The simplest series 
of chain diagrams is 



where 
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and corresponds to Hubbard-I approximation for single-electron Green function. The expres- 
sion 

= G kM 



g l (uj n )-tt ' 



(20) 



in the momentum-frequency representation corresponds to the sum of graphs ([Tj|). The poles 
of function G^(u} n ) determine the spectrum of the single-electron excitations 



e n>l m) = \/2{2E + t % ) =F 1/1 J g 2 + 4t n (S*)g + t 



(21) 



Let us now return to the problem of summation of the diagram series for average (Sf ) taking 
into account the above mentioned arguments. The diagram series has the form 



(Sf) = • 




(22) 



The analytical expressions for loop has the next form 




2 ^ t 



— Y — 

N f g 1 (uj n ) - t% \iuj n - £ ioJ n - £ 

n,k K x 



P7 p; 



= (3( ai P+ + a 2 Pr 
where we used the following notations 



(23) 



2 t\ 1 

2 v tl l__ 

a2 N/3 (<T X K) ~ %) (iw n - e) ' 

Using decomposition into simple fractions and summation over frequency we obtained 

«i = JjY,k( A M £ i(k)) + B x n{e„{t % ))) , 

it 



«2 = ^J2k( A M^(k)) + B 2 n{e n {t % ))) 

k 



where 



At 



eiit-. 



£i(k) - £ n(k) 



B 1 



£n(k) - j 
£n(k) - e i(k) 
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A, 



£/(*r) - £ 



The equation for (Sf) can be presented in the form 



Bo 



£n(h) - si{t- k 



(Sf) = (St) - (St(3{ ai P+ + a 2 Pf)) oc + 



Let us introduce 



Hmf = ^2 



MF 
i i 



where 

H t MF = H io + a x P? + a 2 Pr . 
Then the analytical equation for (Sf) can be expresed in the form 

Sp{Sfe-P HMF ) 



(Sf) = (Sf) 



MF 



Sp(e 



if/? 1 + e _/fe 

- tanh | -(h + a 2 - a x ) + In — 



(24) 



Difference ct2 — «i corresponds to the internal effective self-consistent field acting on pseudospin 

2 ^ £ — £ 



«2 — 



(25) 



4 Mean value of particle number. 

The diagram series for average (rii) (using the perturbation theory, Wick's theorem and expan- 
sion in semi-invariants) can be presented in the form 



(rii) =■=<©- 






(26) 





where 



(P a ) 



lUJr, — £ L 
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(P+;P- 



and the last term appears due to the pairing of the electron creation (annihilation) operators 
with the operator of particle number. 

Analytical expression for (B^) can be obtained starting from formulas flUf) , (|TT)|) 

- <*>» + w s w£)-t ( Hi£U> ■ (27) 

n,/c,a 

where 

Sp{n i e- I3HMF ) 
(n t ) MF = Sp{e - PHMF) ■ 

After simple transformation we obtain next relation 

(n) MF - 71(e) - n(e) = 2(S z )(n(e) - n(e)) , (28) 

or 

(n) M F = 2(P+)n{e) + 2<P~)n(e) , 

where nfe) = — ^-3- is Fermi distribution. 

w l + e^ 

Using decomposition into simple fractions, summation over frequency and relation fl2"8|) we 
can present average (n) in the form 

(»> = I EM £ K%)) + ™M%))} " ^P + )<e) - 2(p-)n(e) . (29) 

k 

5 Thermo dynamical potential. 

In order to calculate the thermodynamical potential let us introduce the parameter A G [0, 1] 
in the initial Hamiltonian 

H\ = H + XH int , (30) 

such that H — > H a for A = and H — > H Q + Hi nt for A = 1. 
Hence 

Z x = Sp(e-^) = Sp{e- pH °a x {(3)) = Z {a x ((3)) , 



where 



and 



r 

a x ((3) = T r exp <-A / H int (r)dT 



n x = -^\nZ -^\n(a x (P)) , (31) 
AQ X = Q X -Q = ~ln(a x ((3)) . 
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Here Q a is the thermodynamical potential calculated with the single-site (diagonal) part of the 
initial Hamiltonian. 
Therefore 



AQ= dX 



o 



(32) 



For the value dQ\/dX, we can write immediately the diagram series in the next form 



P 



dQx 
~~dX 






(33) 



where 



x t and also 



The expression (|32 ) can be presented in the form (using the diagram series (]33f)) 



n,fc 

2 



Xt%g\{u n ) 



N{3 



Eln(l-^K))-— 

n,fc n,fc 



dA. 



The first term in expression ( |34| ) may be written in the diagram form as 



(34) 




(35) 



The series (|35| ) describes an electron gas which energy spectrum is defined by the total pseu- 
dospins field. This series is in conformity with the so-called one loop-approximation. 
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The second term in expression ([53) can be integrated to the following diagram series 




and appears due to the presence of pseudospin subsystem. 

Finally, the diagram series for (3AQ may be written as the sum of expressions (53) and (|36] 
the corresponding analytical expression is the following 



AQ = — — V In (coshf gJ fe))(coshf gjjfe ))_ 
N/3 j (cosh fe) (cosh fe) 



1 f /3 1 + e _/3e 1 

--In cosh | - (ft + a 2 - ai) + In — ^ j + 

1 f /3 1 + e _/3e 1 

+ - In cosh | |/i + In - - e _ p . | + (S z ) (a 2 - ai) . 

Here, the decomposition in simple fractions and summation over frequency were done. 

Then, since the thermodynamical potential is the function of argument (S z ), let us check 
the consistency of approximations made for (S z ), (n) and thermodynamical potential Q. For 
this purpose one should obtain average (S z ) and average (n) from the expression for grand 
thermodynamical potential 

<,{J | £{™M%)) + n{e n {t % ))} - 2<P>(s) - 2(P>( e ) , 



d(-/x) iV 



A- 



- tanh < —{h + a2 — ocx) + In ■ 



d(-h) 2 \2 V ' v 1 + e-^ 
We thus obtain 

f/(> <»>, 3^7 = 



d(—fi) 1 ' d(—h) 

Therefore, the calculation of the mean values of the pseudospin and particle number opera- 
tors as well as the thermodynamical potential is performed in the same approximation which 
corresponds to the mean field one. 

6 Pseudospin, electron, and mixed correlators. 

In this section our aim is to calculate correlators 

K-(T-T') = (TSt(T)S Z m (T')) c , 

KZ(T-T') = (TS!(r)n m (r'))c, 



10 



KZ{r-r') = {Th l {r)n m {r')) c 

constructed of the operators given in the Heisenberg representation with imaginary time argu- 
ment. 

Let us present the diagram series for correlation function (in the momentum-frequency 
representation) within the generalized random phase approximation. In our case (absence of 
the Hubbard correlation) this approximation is reduced, because the so-called ladder diagrams 
(see. [10]) with antiparallel lines disappear. This reduce allow to take into account mean values 
of pseudospin found self-consistently within the mean field approximation. 

We would like to remind that we have neglected diagrams which include semi-invariants 
of the higher than first order in the loop and also connection of two loops by more than one 
semi-invariant. 




where, we define 




P« = (P+ « =(0,1); e a = (e,e); © = P a , O = S* . 

Here, the first term in equation (|38D takes into account the direct influence on pseudospins of 
the internal effective self-consistent field and is given by 




(39) 



Series (|39|) means second-order semi-invariant which renormalized due to the 'single-tale' parts, 
and thus is calculate by H MF . 

The second term in equation ([$8|) describes the interaction between pseudospins which is 
mediated by electrons (energy of electron spectrum is defined by the total pseudospin field). 
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We introduce the shortened notations 



n 



a,(3 





(40) 



Solution of equation ([38|) can be written in the analytical form 

1/4 - (^) 2 



{S z S z )i 



1 + E( _l ) ^^(i_ ( ^ )2) 

»,/3 



where 



n 



k 



T a (k, u ri 



1 



1 



(iu n - e a ) (1 - t n g(u n )) 



(41) 



(42) 



(43) 



Decomposition of the function T a (k, u n ) into simple fractions and subsequent evaluation of the 
sum over frequency leads to the next expression 



a,/3 



* N ^ [ej(t s ) - £//(%)] - s n (t u<r )} 



x 



n[si{t % )\ - n[e n (t % -)] n[e n {t % )\ - 7i[ej(t^)] 1 



k+q) ^Hy^k) ^n^k+qJ ) 

After the substitution ([44]) in equation fl4i"]) we obtain, finally, expression for (S Z S Z ). 
This formula for the uniform case (g = 0) can be rewritten as 

(S*S*) iM > = (l/4-(S*)*)x (45) 

+ 2N f [ ej (t s ) - euik)? { cosh 2 + cosh2 j) l 4 < ; j 

Expression ( [45]) can be obtained from the derivative d(S z ) / d(/3h). This means that mean values 
of pseudospin and pseudospin correlators are derived in the same approximation. 
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For mixed correlator the diagram series has the form 




where 



/= CQ_3><o3>- £ 




and 



n= CO O 




<A + 




=CZD»=«±r» , =Q») 



\3^S— + 



P a T a (k,co r , 




Solution of equation fl47l) can be written in the analytical form 

I = 2(n(e)-n(e))(S z S%. 
Here we start from formula ([y]) and from the next relation 

{S*n) MF -(S*)(n) 



I " <^> 2 



2(n(e) - n(e)) . 



The second term in the diagram series (|4q) we can present as 



a- 



x 



x 



n[e/(*jf)] - n[ £/ (^ +J 0] n[ej(tg)] - n[£//(% +f )] 



ezr (%) - ei (% +(f ) e„ (%) - e n (t 



k+q) 
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Let us introduce the shortened notations for the expression ([51]) 



II=(S z S z ) ? x (52) 

In this way we obtain 

(S z n)t= 2{n{e) - n(e))(S z S z ) (f + (S z S*)fX (53) 

From our diagram series we can see: the correlators containing pseudospin variable S z are 
different from zero only in the static case. This is due to the fact that the operator S z commutes 
with Hamiltonian being the integral of motion. 

For electron correlator our diagram series has the form 




and only last term is not equal to zero for non-zero frequencies. Let us consider series ([54]) 
term by term 



a, P 




(55) 



After simple transformation we can obtain the next relation 



(nn) M F ~ (n) 2 - - 



(P+) 



+ 



(P~) 



2Vcosh 2 ^ cosh 2 f. 



((nS z ) MF -(n)(S z )y 



(56) 



<p+)<p-> 

This relation makes possible to write immediately the simple analytic expression for series (]55|) 



[2(n(s)-n(e))Y(S z S z U+- 



<P+> 



2Vcosh 2 ^ cosh 2 f, 



>5(u) 



Analytical expressions for //-term can be obtained starting from formulae ([48) - (|52|) 

// = {2[n{e) - n(e)]{S z S z )<i x [®]*KM ■ 
Using the expression (|S3) we can unite fl5"5|) and ( |57| ) 



/ + // = { 2 (n(e)-n(e)){S z n) ? +- 



(P + ) 



(P-) 



2 I cosh 2 f cosh 2 f, 



>6{u) 



(57) 



(5f 



(59) 
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The diagram series for the fourth term in (|54| ) has form 




Last term can be presented in the form 






(62) 



Let us the take down final formula for electron correlator for the uniform (q = 0) and static 
(lj = 0) case 

1 ( (P+) (P-) 



(nn) = 2(n(e) - n(i))(S z n) ?=0 + - , ^ ^ 



+ 



COSh 2 f , 



+ 



2JV ^ e/ (t s ) - ejr (t s ) 



1 



cosh 



2 /9ezr(%) 



cosh 



2 /3£/(*e) 



(S"n)^o+ 



l 



l 



2 

2 



cosh 



2 /M*e) 



+ 



2 I cosh 2 ^ cosh 2 f 



(63) 



fc ^ won 2 2 ) v 2 — " 2 ' 

Same result we can obtain from the derivative d(n) /(d/3/i). Thus all our quantities: mean 
values of the pseudospin and particle number operators, thermodynamical potential as well as 
correlation functions are derived in the framework of one approximation which corresponds to 
the mean field approximation. 
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In the second part of the paper we shell perform numerical calculations for the analytical 
expressions obtained in the first part. We shall investigate values of pseudospin and particle 
number operators with the change of the asymmetry parameter h (T = const) or with the 
change of temperature T (h — const) for the cases of the fixed chemical potential value (regime 
H = const) and constant mean value particle number. 
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